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Abstract 

We present results here from 2.5-D particle-in-cell simulations showing that the electrostatic 
(ES) components of broadband extremely low frequency (BBELF) waves could possibly be 
generated by cross-field plasma instabilities driven by the relative drifts between the heavy and 
light ion species in the electromagnetic (EM) Alfvenic component of the BBELF waves in a 
multi-ion plasma. The ES components consist of ion cyclotron as well as lower hybrid modes. 
We also demonstrate that the ES wave generation is directly involved in the transverse 
acceleration of ions (TAI) as commonly measured with the BBELF wave events. The heating is 
affected by ion cyclotron resonance in the cyclotron modes and Landau resonance in the lower 
hybrid waves. In the simulation we drive the plasma by the transverse electric field, E y , of the 
EM waves; the frequency of E y , tOd, is varied from a frequency below the heavy ion cyclotron 
frequency, Qh, to below the light ion cyclotron frequency, Q,. We have also performed 
simulations for E y having a continuous spectrum given by a power law, namely, |E y | ~cc>d’ a , where 
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the exponent a = _, 1, and 2 in three different simulations. The driving electric field generates 
polarization and ExB drifts of the ions and electrons. When the interspecies relative drifts are 
sufficiently large, they drive electrostatic waves, which cause perpendicular heating of both light 
and heavy ions. The transverse ion heating found here is discussed in relation to observations 
from Cluster, FAST and Freja. 

1, Introduction 

In recent years, two quite different low-frequency (FF) wave phenomena have been 
identified in the auroral plasma and they are the electromagnetic ion cyclotron (EMIC) waves 
and the broadband EFF (BBEFF) waves [Boehm et ah, 1990; Bonnell et ah, 1996; Kintner et al, 
1996; Wahlund et al., 1998; 2003; Klatt et al., 2005]. The EMIC waves occur in certain 
frequency bands in a multi-ion plasma [e.g. see Erlandson et al., 1994]. On the other hand, 
BBEFF waves appear without any distinctly identifiable features in the frequency spectrum 
extending from the lowest frequency Alfvenic waves to proton plasma and/or lower hybrid (EH) 
frequency. The low-frequency components in the BBEFF waves are electromagnetic, especially 
below the heavy ion cyclotron frequency in a multi-ion plasma, and they become increasingly 
electrostatic with increasing frequency. For the purpose of easy reference, we call the EM and 
ES components of the BBEFF waves EMEFF and ESEFF waves. For our purposes in this paper 
the EMEFF waves include EMIC waves as well. 

The data from the Freja satellite clearly demonstrated that when EMEFF waves dissipate at 
altitudes near h~ 1700 km, some of the wave energy goes in transverse heating of ions [Fouarn 
et al., 1994; Wahlund et al., 1994; Erlandson et al., 1994], Erlandson et al. [1994] reported 
interesting examples of EMIC waves in two bands, one below the 0 + cyclotron frequency and 
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the other between He + and H + cyclotron frequencies, the latter correlating with preferential 
heating of 0 + ions. Results from FAST also show strong connection between EMELF waves 
and ion heating [e g., see Chaston et al., 2004]; the ions’ pitch angles measured from FAST at 
altitudes near h~3500 km indicate local transverse heating and most energetic ions are the heavy 
0 + ions. 

In a series of recent satellite observations from the Polar satellite, it is demonstrated that 
plasma electrodynamics along the field lines threading through the pole -ward boundary of the 
auroral oval and the plasma sheet boundary layer is dominated by Alfven waves [ Wygant et al., 
2000, 2002; Keiling et al., 2005; Angelopoulos et al., 2002]. Specifically Wygant et al. [2002] 
show that narrow-structured Alfven waves, believed to be the kinetic Alfven Waves (KAWs) are 
embedded in much wider Alfven wave structures, which have long wavelengths across the 
ambient magnetic field. It is reported that such waves accelerate electrons parallel to the ambient 
magnetic field and ions perpendicular to it. 

Wahlund et al. [2003] have reported observation BBELF waves from Cluster at a high 
geocentric altitude of ~ 4-5 Re showing that features of the wave are similar to those seen from 
rocket and low-altitude satellites [Boehm et al., 1990; Wahlund et al., 1994; 1998; Stasiewicz et 
al. 2000]. Figure la shows the frequency spectrum of the fluctuations in the perpendicular 
electric field of the ELF waves. The two curves in the figure, one with lower power level in a 
quiet region and the other curve in an immediately following active region with intense waves 
during the course of traversal of Cluster show the increase in the electric field wave power by 
one to three orders of magnitude in the active region over the entire band of frequency shown. 
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Wahlund et al. [2003] point out that for the spatial scales of the waves, the Doppler shift due to 
satellite velocity is not significant, unlike in the case of Freja [Stasiewicz et al., 2000] and FAST 
[Chaston et al., 2004] observations. So the frequencies shown are in the rest frame of the plasma. 
The power in the magnetic component during the intense wave period increases below 2 Hz, and 
only below 0.4 Hz as much as the electric power (Figure lb). Although the spectrum in Fig. 2a is 
only up to 10 Hz, Wahlund et al. [2003] point out that the enhancement extends up to the low- 
frequency cut-off of the auroral hiss emissions observed near 300-400 Hz. Furthermore, they 
found that 6E/5B steadily increases above 0.4 Hz, and below this frequency it is close to 9xl0 6 - 
3xl0 7 m/s, the range of the Alfven velocity during the observation time. It is worth mentioning 
that 0 + cyclotron frequency is 0.46 Hz and for H + it is ~ 7.5 Hz. Thus, the low-frequency 
components below ~0.4 Hz are Alfvenic waves while at higher frequencies the waves become 
increasingly electrostatic. Klatt et al. [2005] recently reported observations of ELF waves, 
showing a direct relation between the Alfvenic waves and electrostatic noise; the latter coincides 
with the peak in the transverse electric and magnetic fields of the former. Since the peak of the 
magnetic field perturbations coincides with the maximum gradient in the current density 
associated with the DAWS, Klatt et al. [2005] have suggested that the noise, which extends both 
below and above the lower hybrid frequency, might be generated by current- advective shear- 
driven CASDI instability (CASDI) [Seyler and Xu, 2003]. Even the Observation from Freja 
shows a direct connection between EMELF waves and high frequency electrostatic waves 
[Wahlund et al., 1994; Erlandson et al., 1994]. Wahlund et al. [1998] suggested that as the 
EMELF waves propagate down, they emit the electrostatic waves. They also suggested that the 


electrostatic waves could be ion-acoustic waves. 
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Khazanov et al. [1997a,b] have theoretically demonstrated that the Alfven wave-driven drifts 
of electrons and ions in a multi-ion plasma could generate high-frequency electrostatic waves 
even when electron and ion temperatures do not satisfy the condition for the existence of ion- 
acoustic waves, namely, T e >3Ti. The high-frequency mode generated by this process is lower 
hybrid waves, which heats ions transversely and electrons parallel to the ambient magnetic field. 
Singh and Khazanov [2004] showed by means of 2.5-D PIC simulations, that Alfven wave- 
driven drifts could generate broadband electrostatic turbulence covering range of frequency from 
the ion cyclotron harmonics to the LH frequency depending on the relative densities of heavy 
and light ions. For low densities of the heavy ions, the electrostatic turbulence develops a distinct 
peak near the LH frequency. As the heavy ion density increases, the peak broadens covering the 
ion cyclotron hannonics. 

The purpose of this paper is to extend the work of Khazanov et al. [1997a,b] and Singh and 
Khazanov [2004] by considering a range of frequencies of the EMELF waves driving the 
plasma. First we consider monochromatic coherent waves, whose frequencies vary from below 
the heavy ion cyclotron frequency £2h to below the light ion cyclotron frequency Q, in different 
simulations. In view of the broad broadband EMELF waves seen in satellite observations 
[Erlandson et al., 1994; Wahlund et al., 2003], we also consider driving fields having a power 
law spectrum. In the latter case we find that the impulsive nature of the driving electric field in 
the time-domain, generates large amplitude electrostatic waves causing a large transverse ion 
heating. 


2. Method of Simulation 
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We use a standard 2.5-D particle-in-cell code [Birdsall and Langdon, 1985], which we have 
described in previous studies [Singh and Leung, 1999, Singh, 2000]. Specifically we simulate 2- 
dimensional plasma of dimension L x xL y in the x-y plane. The ambient magnetic field B 0 has only 
the z- component in all the simulations reported here. In a set of three simulations from R-l to 
R-3, plasma is driven by a monochromatic sinusoidal electric field given by 

E y = a y F a (t) EoSin(wdt), F a (t) =erfc( t/x s ) , (1) 

in which the driving wave frequency oc>d is varied from run to run (R1-R3) ranging from below 
the heavy-ion cyclotron frequency £2h to just below the light-ion cyclotron frequency Q;. The 
plasma has a heavy ion species along with the light ions; the relative density of the heavy ions is 
10%. The driving field amplitude E 0 is mentioned while discussing the simulation runs. F a (t) is 
an adiabatic start function which allows E y (t) to grow slowly over a time period x s to reach the 
steady state amplitude E 0 . x s is several times the driving wave period of the heavy ions. If x s is 
not sufficiently long or when F a (t)=l, the response of the ions contain spurious oscillations. In 
order to economize computing time, we have used ion masses which are not real; we have used 
Mi/m=400, where M; is the light ion mass and m is electron mass. For the heavy ion mass, we 
used Mh/Mi=4. In the other set of three simulations (R4-R6) the monochromatic driving electric 
field is replaced by a broadband spectrum covering the frequency range from below £2h to ff . 
The features of the power-law frequency spectra in these simulations are discussed later. The 
size of the simulation box is L x xL y = 256x256Xd 2 . We employed 25 electron-ion pairs per cell 
with cell size XdxXd, and Q e /o)p C = 1 , where ) H \ is the Debye length, and Q e and o) pe are the electron 
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cyclotron and plasma frequencies, respectively. We used a time step At=2jt/32. In the PIC code 
and our following discussions we use nonnalized quantities as follows: distances are normalized 
in the units of Xj, velocity in units of electron thermal velocity a eo , time in units of the inverse of 
electron plasma frequency co pe o, electric potential in units of cp n =T 0 /e, and electric field in units 
of En = _ n |_d, where X ( i=a c /o) pc , and a eo =(T 0 /m) . Ion thermal velocities are likewise defined by 
a 0 =(T /M ) , where _ denotes the ion species. T 0 is the initial plasma temperature, and it is the 
same for the electrons as well as the light and heavy ions. 

3. Numerical Results 

First we present results from simulations, in which the driving EMELF waves are 
monochromatic. Since the EMELF components of EMIC and BBELF waves overlap in 
frequency spectrum, such simulations could be relevant to TAI by both EMIC as well as BBELF 
waves. The monochromatic simulations are followed by discussion of results from simulations 
driven by broadband EMELF waves. We not only discuss here TAI, but also the generation of 
ES waves, which are an integral part of the events of EMIC waves and of course BBELF waves 
observed from satellites. 


3,1 Monochromatic Driving Wave 

The results from run R-l are physically analogous to previously reported simulations [Singh and 
Khazanov, 2004], except that heavy ion mass is Mh=4M;, in contrast to Mh=9M, in the previous 
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study. In R-l we have wave frequency co= 0.75Qh=0. 1 875Q,. Figure 2a shows the temporal 
evolution of the wave amplitude as seen in the electrostatic potential cp (t) at the mid point (L x /2, 
Ly/2) of the simulation region. The unit of cp is T e /e and the time on the horizontal axis is labeled 
in units of light ion cyclotron period x c i. The evolution of perpendicular temperature of the heavy 
ions, Tj„ and of the light ions, T ± i, as well as of the parallel temperatures of the ions and 
electrons are plotted in Figure 2b. Since the ambient magnetic field is purely perpendicular to the 
plane of the simulation, as expected the parallel temperatures remain constant at the nonnalized 
value of unity. As the driving field E y (t) grows in the simulation according to the scheme of 
adiabatic start (see eq. (1)), the amplitude of noise in the electrostatic potential grows. This 
growth is caused by the instabilities generated by the relative drifts between light and heavy ions 
in the driving field E y . As an example, the evolution of the ExB drift, V x h, and the polarization 
drift, V y h, of the heavy ions are plotted in Figure 2c. When E y (t) reaches its steady oscillation as 
E 0 sin((nt), the waveforms of the drift velocities match the theoretical drifts given by 

V xs =[l/(l-co 2 /Q s 2 )] (Eo/Bo) sin(o)t), (2) 

V ys =[(to/Q s ) /( 1 -or/a 2 )] ((Eo/B 0 ) cos(cot), (3) 

where the subscript ‘s’ stands for the ion species with s=h(i) for heavy (light) ions. In the steady 
state, the amplitudes of V x h and V y h are 0.57a eo and 0.42 a eo , respectively, and they match well 
with those given by the above relations. The light ion drift amplitudes are V XI ~ 0.25 a eo and 
V y i=0.047 a eo . For the electrons, the amplitudes are V xe = 0.25 a eo and V ye ~ 0. As found before, 
mainly the heavy and light ion coupling due to the relative polarization drift, V yre i= V y h-V y i = 0.2 
a eo = 4a l0 , drives the instability to generate the fluctuation in cp(t) plotted in Figure 2a. 
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Figure 2b shows that until t=23x c ;, T ± h and T ± ; increase more or less smoothly, then there 
are two sudden jumps in T ± h, separated by one driving field wave period x 0j . The step increases in 
the temperature coincide with bursts of large amplitude lower hybrid waves (Figure 2a). The 
frequency spectrum of the data plotted in Figure 2a is shown in Figure 2d. The noteworthy 
feature of the spectrum is the peak appearing near the lower hybrid frequency o>ih -13.6Q,; the 
lower hybrid frequency is given by 

com = tOpifn; +(n h /4)] 1/2 x (l+oo pe 2 / Q e 2 )Y m , (4) 

where to p i is the light ion plasma frequency with the total density n 0 =nh+ni and it is o) pi ~20Q 1 for 
Mi/m=400. Since Q e /o)p C = 1 , nih = 0.1n o and ni=0.9n o , we have coih= 1 3.6Q,. We point out that 
even though the spectrum has a peak at coih, it is quite broadband extending even below the heavy 
ion cyclotron frequency. The low-frequency behavior of the spectrum is perhaps an artifact of 
the adiabatic start in which the driving wave grows over several ion cyclotron periods. But the 
reason for wave power between the heavy ion cyclotron frequency and ooih is the generation of 
waves by the drifts shown in Figure 2c. 

We point out that the LH waves are driven by the interaction of the slow-beam mode of the 
heavy ions and the Bernstein mode of the light ions [Seiler at ah, 1976]. For such interactions the 
ion-ion instability growth rate maximizes near the light ion cyclotron harmonic just below the 
lower hybrid frequency, i.e., oiih-n Qj; in R-l discussed here n=13. Seiler et al. [1976] studied 
this ion-ion interaction by linear instability analysis as well by laboratory experiment. The same 
ion-ion interaction mechanism operates when pick-up ions interact with the background plasma 
generating LH waves [Cairns and Gumett, 1990; Singh, 2000]. Both heavy and light ions diffuse 
in velocity space under the influence of this instability affecting their heating. It is worth noting 
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that for such ‘high-frequency’ ion-ion interactions the heavy ions behave as un-magnetized ions 
because the interaction time is much shorter than the heavy ion cyclotron period. When the LH 
waves dominate the heating occurs primary by Landau resonance co-kyV,, where k y is the 
perpendicular wave number and V; is the ion velocity. 

We perfonned runs R-2 and R-3 with driving wave frequencies to/ £2, =0.5, and 0.8, 
respectively. Figure 3 shows the comparison of the ion heating in the three runs R-l to R-3. The 
top panels show the V x -V y phase space of light (black dots) and heavy (red dots) ions at initial 
time t=0. The middle panels show the same phase space at the end of the simulations (t=65.62 
x c i). A noteworthy feature of the phase space in the middle panels is that they are not centered at 
V x =V y =0; they are off-centered due to the combined effects ofExB and polarization drifts. The 
bottom panels show the evolution of the heavy and light ions’ perpendicular temperatures, like in 
Figure 2b. The frequencies of the driving waves are labeled in the top panels. It is interesting to 
note that for the same amplitude, the driving wave becomes more effective in heating the ions as 
the wave frequency moves closer to the light ion cyclotron frequency. It is seen from the much 
larger spread of the phase space for co/ £2; =0.8 than that for co/ £2; =0. 1875 and 0.5. We also see 
that heating rates for co/ £2i =0.8 are much larger than that for co/ £2; =0. 1875 and 0.5. 

The reason for the relatively large heating rate as co tends toward the light ion cyclotron is the 
reduction in the drifts of heavy ions and the increase in the light ions drifts due to the resonance 
effects detennined by the factor (l-co 2 /£2 s 2 )' 1 in (2) and (3). Figure 4 shows the evolution of the 
drift of both light and heavy ions as well as of electrons from R-3, in which co=0.8£2i. The solid 
and dotted lines show the ExB and the polarization drifts respectively. The ExB drift for the 
heavy ions is V x h = 0.024a eo , while for the light ions it is V x i= 0.7 a eo giving a relative drift 
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V x , re i= 0.72a eo . Note that the individual ExB drifts are out of phase, so they add to yield the 
relative velocity. The polarization drifts are V y h=0.09 a eo and V y i=0.55 a eo , resulting into a 
relative drift V yjre i=0.63 a eo . Note that V x . re i is slightly larger than V y>re i- 

Now let us consider the drifts if we had taken co=0.8£2h, below the heavy ion cyclotron 
frequency, the in-phase ExB drift of heavy and light ions are V x h= 0.7 a eo and V x i= 0.25 a eo , 
respectively, yielding a relative drift V x , re i2~0.45 a eo , which is much smaller than the above 
estimate of V x , re i. The polarization drift of heavy and light ions are V y h= 0.55 a eo and V y i= 0.05 
a eo , yielding a relative drift V y , re i 2 =0.5 a eo , which is also smaller than the V y> rc i calculated above. 
Thus we learn that the resonance effect, in combination with the drifts of heavy and light ions 
being out-of-phase as the driving wave frequency increases above the heavy ion cyclotron 
frequency, produces much larger relative drifts between the light and heavy ion species than that 
for frequency being below the heavy ion cyclotron frequency. The larger the relative drift, the 
larger becomes the amplitude of the driven electrostatic waves and the resulting ion heating. This 
might explain the observation from Freja that the transverse ion heating correlates better with the 
EMIC waves above the He + cyclotron frequency than that with the waves below the 0 + cyclotron 
frequency as reported by Erlandson et al. [1994], 

Figure 5a shows the evolution of the noise in the electric potential from R-3 as in Figure 2a 
for R-l. The wave amplitude in Figure 5a is generally larger than that in Figure 2a. Figure 2b 
shows the temporal behavior of the waves on an expanded time scale revealing that they contain 
strong features at light ion cyclotron period and the periods corresponding to higher cyclotron 
hannonics. The frequency spectrum of the noise in Figure 5a is plotted in Figure 5c. Comparing 
the spectra in Figure 2d and Figure 5c, we find that in the latter case the peak at the lower hybrid 
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(LH) frequency is washed out and the spectrum is nearly flat from a low frequency near the 
heavy ion cyclotron frequency (0.25£2i) to the LH frequency, o)ih~13.6Q,. Furthermore, the 
power in the flat part of the spectrum in Figure 5 c is about one order of magnitude larger than 
that in Figure 2d. A detailed examination of the wavefonns in Figure 5a reveals that the noise 
consists of ion cyclotron oscillations at light and heavy ions’ cyclotron frequencies as well as 
high frequencies extending up to the LH frequency. We leam from this simulation that EMELF 
wave electric fields near the light ion cyclotron frequency belonging to EMIC waves (Erlandson 
et al, 1994] as well as BBELF [Wahlund et ah, 1998] waves could potentially generate 
broadband electrostatic turbulence. 

3.2. Broadband Spectrum of the Driving Wave 

In the above discussion we have examined the heating by monochromatic driving waves. 
However, observations reveal that BBELF waves are generally broadband extending in 
frequency from below the heaviest ion cyclotron frequency to above the light ion cyclotron 
frequency. Such waves at different frequencies might not be correlated either in time or space. 
The amplitude of such waves generally decreases with the increasing frequency, with exceptions 
of some peaks in the power spectrum. We model such waves as superposition of a large number 
of frequency components with random phase and amplitudes as follows: 


E y (t) = F a (t) G(t), 

(5a) 

G(t) = Z N j=i Aj sin(o)jt +0j), 

(5b) 

Aj= Ai (ooi/coj) a , 0j=Rj 2 7t, 

(5c) 


where the exponent a=l in R-4, a=2 in R-5 and a=l/2 in R-6. The phase 0j is uniformly 
distributed between zero and 2 it detennined by the random number Rj uniformly distributed over 
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0<Rj<l. F a (t) is the adiabatic start function defined earlier. N is the number of mono-chromatic 
waves in the broadband spectrum. We have used Ai=0.25, and o)i=0.1£2j. The separation 
between the consecutive frequencies is iOj+i- ooj =0.01 IQ, with j ranging from one to N=100. 

Thus the amplitude of the component with the maximum frequency is Ai/10=0.025 in Figure 6a 
(R-4) and Ai/100=0.0025 in Figure 8a (R=5). Figures 6a and 8a show G(t) for a=l and a=2, 
respectively. We notice from Figures 6a and 8a that the superposition of several sinusoids with 
random phase produces impulsive electric field behavior. The impulses have large amplitudes 
several times the largest amplitude Ai. The smaller the exponent a the narrower and larger the 
impulses become. For _ =1/2, the impulses become unrealistically large and we do not discuss 
any results from R-6. How do the impulses affect the generation of electrostatic wave and 
associated ion heating? 

Figure 6b shows the time history of the electric potential at the point (64, 64), like in Figure 
2a from R-l. Figure 6c shows the corresponding evolution of the heavy and light ions 
perpendicular temperatures in the same format as in Figure 2b. Note that coinciding with the 
burst of electric noise (Figure 6b) near t~6-7x c i the perpendicular temperatures Tj, and T ± ; rise 
sharply; the heavy ions’ temperature T ± h rises to ~45T 0 . while for the light ions the increase is to 
Tj_i =15T 0 . After the sudden jumps in the temperatures the heating is relatively slow. The light 
ions heat with a nearly constant heating rate. The heavy ions heat discontinuously over a certain 
time range, but the jumps in the temperatures are relatively small. The sharp rise in temperatures 
near t~6-7x c i coincides with an impulse in the same time frame marked in Figure 6a by fi. The 
small jumps in the heavy ion temperature in the time frame t~ I 4-20 t ci are associated with weaker 
bursts of noise seen in Figure 6b. When the ions are heated, the noise bursts becomes relatively 
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weak too as expected from theoretical considerations; the threshold criterion on the relative drifts 
are not met after the heating of the ions [Singh and Khazanov, 2004]. It is noteworthy that the 
large-amplitude impulses in E y (t) in the time frame from t > 30x c , (Figure 6a) do not affect 
further heating of the heavy ions while the heating of the light ions continues. This is attributed 
to the fact that preferentially heated heavy ions have large Larmor radii and the ES waves 
heating them have comparable wavelengths. Thus the heating saturates for the heavy ions, while 
for the light ions with smaller Larmor radii the heating continues slowly till the end of the 
simulation. 

Figure 6d shows the frequency spectrum of the electric potential plotted in Figure 6b. Note 
that the electrical noise becomes broadband extending up to the lower hybrid frequency 
coih= 1 3.6Q,. The spectrum contains peaks near the heavy ion cyclotron frequency and its 
hannonics. Comparing the spectra in Figure 2d from R-l and in Figure 6d, we find that the 
amplitude level at the frequencies below 10 (On, is enhanced by an order of magnitude in 
R-4. Furthennore, the spectrum has become flatter. 

In order to understand the origin of the large amplitude burst near t~6-7x a , we examine the 
drift velocities of the light and heavy ions in this time frame. Figures 7a-7b show V x (ExB drift 
Vd) and V y (polarization drift Vd) for the light and heavy ions, respectively. The differences (V x i- 
V xh ) and (Vyi-Vyh), giving the relative drifts in the x and y directions between light and heavy 
ions, are plotted in Figure 7c and 7d, respectively. Note the largest relative drifts occur for the 
first time in the period t~6-7x c i. The drifts coincide with the impulse f in Figure 6a; the relative 
drifts caused by ExB drift and polarization drifts are nearly of equal magnitude ~ 0.8 ae 0 ~ 16 ai 0 . 
Such drifts are large, but occur over a duration shorter than the ion cyclotron period, and produce 
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large-amplitude bursts of electrostatic noise. At latter times even larger drifts appear (not 
shown), but they do not generate similar large amplitude burst as in the time period t~6-7x c i. The 
reason for this is that the ions heated in the earliest event of the large-amplitude burst require an 
even larger drift for further wave generation [Singh and Khazanov, 2004]. 

Like in Figure 6a-6d for R-4, Figure 8 shows the time history of (a) E y (t) for a=2 as 
already mentioned, (b) electric potential cp(64,64), (c) ions’ perpendicular temperatures and (d) 
the frequency spectrum of cp(64,64) from R-5. The frequency components of the driving field 
E y (t) in R-5 decay faster with the increasing frequency than that in R-4. This results into smaller 
amplitude of the impulses in E y (t) as seen in Figure 8a. In turn, this results in smaller amplitude 
bursts and at a later time in R-5 than that in R-4. The ion heating shows the steps in the 
perpendicular heating of both heavy and light ions near t~20r cl (Figure 8c) , but the step sizes in 
Tj, and T , are relatively smaller in R-5 than that in R-4. Each burst in cp(t) is associated with a 
clearly visible step in heating of the heavy ions over the interval t~20-30r cl . The frequency 
spectrum in Figure 8d shows the LH peak and additional peaks below the LH frequency. 

The appearance of the LH peak in Figure 8d and its absence in Figures 6d and in Figure 5c 
have the following physical implication. When a high frequency EMELF wave dominates the 
frequency spectrum of the BBELF waves, the drift generated ES waves tend to have flatter 
spectrum than that for the EMELF spectrum being dominated by the low frequency waves. 

4. Conclusions and Discussions 


Our main conclusions are the following. 
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a. When an EMELF wave permeates a plasma, the drifts of the charged particles in the 
transverse electric field of the wave generate electrostatic waves, which heat the ions 
transversely. 

b. The frequencies of the generated electrostatic waves range over the ion cyclotron harmonics 
to lower hybrid frequency, appearing as broadband electrostatic noise. 

c. When EMELF waves constitute a broad spectrum including electromagnetic ion cyclotron 
waves ranging in frequency from below the heavy ion cyclotron frequency up to the light ion 
cyclotron frequency, their superposition produces impulsive electric fields of large 
amplitudes in time domain. 

d. The drifts in the large-amplitude impulses produce a burst of large-amplitude electrostatic 
noise, which causes sudden perpendicular heating of both heavy and light ions. For the heavy 
ions several events of such step heating occur while the light ions undergo a smooth heating 
after the first step. 

e. For the amplitudes of the waves considered here, the perpendicular heating is significant. The 
ions can be heated up to several tens of the initial ion temperature. 

f. For the EMIC waves, the waves having frequencies between the light and heavy ions are 
more effective in transverse ion heating due the larger relative drifts between the different ion 
species due to the 180 phase difference between the drifts of the individual species. This 
phase difference comes from the resonance factor 1/(1 -co 2 ,/Q 2 ); for heavy ions this factor is 
negative while it is positive for the light ions. 

g. The flatness of the generated ES wave spectrum depends on the dominance of the high 
frequency component of the EMELF (Alfvenic) waves. 
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Our simulations, as discussed here, demonstrate that EMELF, including 
electromagnetic ion cyclotron waves, can be dissipated via the generation of electrostatic waves 
produced by the relative drift between ions of different species by the former waves. A direct 
application of our results presented here is the observation of broadband extremely low 
frequency waves (BBELF) and associated ion heating. BBELF waves have been observed in 
rocket experiments [Boehm et ah, 1990], and from several satellites such as Freja [Fouarn et ah, 
1994], FAST [Stasiewicz et ah, 2000; Chaston et al, 2004] and Cluster [Wahlund et ah, 2003]. 
However, it is now well established that EMEFF waves are invariably associated with high- 
frequency electrostatic waves extending in frequency up to the lower hybrid frequency, which is 
close to the H + plasma frequency. These assertions are highlighted by the Cluster observation of 
BBEFF waves reported by Wahlund et al. [2003] and more recently by Klatt et al. [2005] based 
on SIERRA observations. As presented by Wahlund et al. [2003], the Cluster data reveals the 
following properties of BBEFF. The observations were made at a high altitude of 4-5 Re. The 
frequency spectrum of both EMEFF and electrostatic components are plotted in Figure 1 up to 
10 Hz and as mentioned by Wahlund et al., it continues with enhanced power level up to 300- 
400 Hz, which is about the H + - plasma frequency f p h =370 Hz. The other characteristic 
frequencies in the plasma for the Cluster observations are proton cyclotron f cl ~7.5 Hz, electron 
plasma frequency f pe — 1 5 kHz, and electron cyclotron frequency f cc ~ 1 4 kHz. For these parameters 
the lower hybrid frequency is fih~260 Hz, close to the approximate upper frequency limit of the 
electrostatic component of the BBEFF waves. Klatt et al. [2005] shows that peaks in the burst of 
electrostatic waves extend to the FH frequency and beyond it as VFF waves. They interpret their 
results to the presence of lower-hybrid solitary structures (FHSS) [ FaBelle et al., 1986; Kintner 
et al., 1992; Vago et al., 1992; Fynch et al., 1999; Schuck et al., 2003]. 
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Despite a large number of studies on LHSS, the energy source for the LHSS remains 
undetermined. A myriad of electrodynamic processes consisting of Alfven waves solitary 
structures, transverse ion heating, formation of density cavities, large parallel electric fields in 
the cavities, field-aligned burst of electrons, transversely limited wave packets of lower hybrid 
waves in cavities (LHSS), and all on the same auroral field lines have been observed. It is 
established now that the transverse heating of ions is caused by Alfven waves as well as lower 
hybrid waves [Andre et ah, 1998; Lynch et ah, 1999]. It is also demonstrated by modeling that 
density cavities are a natural consequence of transverse heating of ions [Singh, 1992; Singh and 
Chan, 1993; Singh, 1994,1996]. Connection between density cavities and parallel electric fields 
in the cavities in association with Alfven waves are established by both observations [Chaston et 
ah, 2006] and modeling of current-cavity interaction [Singh, 2002; Singh, 2005], The parallel 
electric fields in the fonn of a double layer in the cavity create an electron beam. The lower 
hybrid waves can be generated by the Alfven waves directly as studied in this paper as well as by 
the electron beam, especially if the beam is thin in the transverse direction [Singh et ah, 1985; 
2001]. However, the knowledge of the electrodynamic processes mentioned here exists only in 
piecemeal. A self-consistent kinetic modeling at meso-scale size is needed to establish the 
mutual coupling between them. Such a study is likely to reveal the source of energy for the 
LHSS events. 

For the EMELF waves at frequencies below the heavy ion cyclotron frequency the 
polarization drift of heavy ions, V y h a given by (3), drives the instability. This drift depends on 
the frequency via the factor H(a) = a/(l- a 2 ) where a= o)/Qh . For a<0.6, the factor H(a) 
reduces the drift rapidly as the frequency decreases. On the other hand, for a>0.6 this factor 
increases rapidly as the frequency increases toward the heavy ion cyclotron frequency. Thus the 
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latter range of the frequency spectrum of the EMELF waves have the potential for driving strong 
instabilities and ion heating depending on the perpendicular electric field because V y h- E /B 0 . 
Since EMELF waves have very small parallel wave numbers they are difficult to have cyclotron 
resonance with thennal ions, especially at low altitudes of Freja and FAST where characteristic 
ion energy is -0.15 eV. The EMELF waves in the frequency range 0.6< 0 )/Qh< I are likely to 
drive the instabilities studied in this paper and heat the cold ions even when they are not directly 
in cyclotron resonance with the waves. Our study shows that if the waves are broadband, the 
resulting time-domain impulses in E ± have large amplitudes, resulting in larger drifts and heating 
of ions (Figures 7 and 8). 

We give some estimates here of the heating based on wave amplitudes and frequencies 
observed from Freja, FAST and Cluster. Since the polarization drift scales as V y h~ H(a) (E_ /B 0 ), 
and instability threshold is V y h~ a;, the light ion evolving thermal velocity [Khazanov et ah, 
1997a], the extent of light ion heating scales as T^-0.5 M, H(oc) 2 (E ± /B 0 ) 2 [Singh and 
Khazanov, 2004], Heavy ions heat somewhat higher. For the Cluster data in Figure 1, E_~ 1 00 
mV/m and B o ~160 nT, giving an estimate of the heating T j~0.5 H(a) 2 M (E_ /B 0 ) 2 -3.5 keV 
with a -0.7. At the altitudes covering Freja and FAST, we have E ± ~ 25-1000 mV/m and 
B 0 ~12, 000-25, 000 nT. Specifically, we take the example from Figure 1 of Andre et al. [1998]; In 
this case B~25,000 nT and E ± - 200-600 mV/m, giving V y h -10-30 km/s for a -0.7. For such 
drifts light ion heating is estimated as 1-10 eV. If we choose a higher frequency, the heating 
estimate is increased; for a -0.9, we have an increase by a factor of -20 giving the estimate -20- 
200 eV. Ion heating within the ballpark of this estimate is reported by Andre et al. [1998]. This 
increase is simply accounted for by the increase in H(a) 2 . Chaston et al. [2004] also report 
temperatures of ions observed with Alfven waves in this ballpark using FAST data. Knudsen and 
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Wahlund [1998] reported micro-bursts of core ions heated to a few eV associated with EMELF 
waves with low amplitudes -25 mV/m. Even for such low amplitudes, the polarization drift of 
heavy ions can be sufficiently high to excite instabilities in H + - O ' plasma and heat the ions; For 
example, if a ~0.9 the drift V y h- 5 km/s reaching the thermal velocity of 0. 1 5-eV H + ions. 

Heavy ions are expected to heat -2.5 eV in such instabilities. As pointed out earlier, if the waves 
are broadband, the resulting time-domain electric field impulses have large amplitudes and affect 
stronger heating compared to the estimates given here for the monochromatic waves. The 
heating by the instability mechanism does not require direct ion cyclotron resonance with the 
EMELF waves. 

The heating by ion cyclotron resonance as proposed by Chang et al. [1986] and applied to 
Freja data runs into several difficulties. First of all Andre et al. [1998] point out that it takes only 
about 1 % of the observed wave power to achieve the measured level of heating. A question 
arises as to why the heating is so severely limited. The other difficulty is determination of the 
power spectral density (PSD) from the measurements. The critical PSD at the 0 + cyclotron 
frequency 25 Hz (-10 4 m 2 V 2 /m 2 /Hz) for ion heating found from Freja is derived from time 
series measured during the traversal of the EMELF waves by the satellite [Knudsen et al., 1998], 
It is not certain whether the satellite measured spatial or temporal features or their combinations 
[Stasiewicz et al., 2000; Chaston et al., 2004], Thus the PSD values estimated from Freja 
measurements need to be reexamined and resulting ion heating reassessed. 

It is important to point out that if ions are already hot having temperatures as estimated 
above, further heating by the instability mechanism is not operative. A direct consequence of this 
finding is that ion heating might occur when the EMELF waves propagate down where ions are 
cooler. At very low altitudes, where electron skin depth and ion Larmor radii become 
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comparable, the EMELF waves might directly heat the ions by demagnetizing them as suggested 
by Cole [1976] and seen in numerical simulations of V-shaped double layers [Borovsky, 1984; 
Greenspan, 1984; Singh et ah, 1987; Singh and Khazanov, 2003] and in heating of ions by 
narrow scale shear Alfven waves [Stasiewicz et ah, 2000; Chaston et ah, 2004]. 

We point out that for a quantitative one-to-one comparison between the features of EMELF 
waves, their continuation to high frequencies as electrostatic waves and associated ion heating as 
seen in satellite and rocket observations and simulations, the latter needs to be performed with 
specific plasma conditions of the observations and real masses of ion species. We also need to 
include in the simulations the finite perpendicular wavelengths of the EMELF waves. We plan to 
do this in the near future. The work in this paper is meant to reveal the basic physical processes 
involved in generation of the electrostatic waves and resulting ion heating when the plasma is 
driven by a single monochromatic EMELF wave or a broad spectrum of them. 
Acknowledgement: This work was supported by NASA grant NAG5 13489 and a grant from 
Marshall Space Flight Center. 
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To ensure prompt publication: 

Figure Captions: 

Figure 1, Observed power spectrum of BBELF waves from Cluster [Wahlund et ah, 2003]: (a) 
Power spectrum of electric field in quiet (bottom curve) and disturbed (top curve) regions 
traversed by Cluster, (b) Magnetic field spectra corresponding to the spectra in (a). 

Figure 2. Temporal evolution of (a) electrostatic wave amplitude, (b) perpendicular ion 
temperatures of light and heavy ions, (c) Polarization (V y h) and ExB (V x h) drift velocities of 
heavy ions, and (d) frequency spectrum of the drift-driven electrostatic waves in (a) from R-l. 

Figure 3. Comparison of transverse ion heating of light and heavy ions as seen from R-l to R-3: 
top panels show the initial V x -V y phase space for the driving frequency labeled in the panels, 
middle panels show the same phase space at the end of the simulations at t=65.62x c i. The bottom 
panels show the evolution of the perpendicular temperatures of heavy (solid line) and light 
(dotted line) ions. 


Figure 4. Temporal evolution of the drifts of light ions (top panel), heavy ions (middle panel) 
and electrons (bottom panel) from R-3, in which (o ( |/Qi = 0.8. 
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Figure 5. Temporal evolution: (a) electrostatic wave amplitude, (b) Same as (a) but over an 
expanded horizontal scale, (c) Frequency spectrum of data plotted in (a) from R-3. 

Figure 6. Temporal evolution of (a) the driving electric field amplitude E y (t) for the power-law 
spectra with the exponent a=-l. (b) the wave amplitude, (c) perpendicular temperatures of 
heavy and light ions and (d) frequency spectrum of the data plotted in (b) from R-4. 

Figure 7. Temporal evolution of the drifts of (a) light ions, (b)heavy ions (c) relative drifts 
between heavy and light ions due to their ExB drifts, (d) same as (c) but due to their polarization 
drifts, and ((e) drifts of electrons from R-4. The time period shown in these panels cover the 
large-amplitude burst of electrostatic noise in Figure 6a. 

Figure 8. Same as Figure 6, but from R-5 in which the exponent in the power law spectrum is 
a=-2. 
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Figure 1. Observed power spectrum of BBELF waves from Cluster [Wahlund et ah, 2003]: (a) 
Power spectrum of electric field in quiet (bottom curve) and disturbed (top curve) regions 
traversed by Cluster, (b) Magnetic field spectra corresponding to the spectra in (a). 
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Figure 2. Temporal evolution of (a) electrostatic wave amplitude, (b) perpendicular ion 
temperatures of light and heavy ions, (c) Polarization (V y h) and ExB (V x h) drift velocities of 
heavy ions, and (d) frequency spectrum of the drift-driven electrostatic waves in (a) from R-l 
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Figure 3. Comparison of transverse ion heating of light and heavy ions as seen from R-l to R-3: 
top panels show the initial V x -V y phase space for the driving frequency labeled in the panels, 
middle panels show the same phase space at the end of the simulations at t=65.62x c i. The bottom 
panels show the evolution of the perpendicular temperatures of heavy (solid line) and light 
(dotted line) ions. 






Figure 4. Temporal evolution of the drifts of light ions (top panel), heavy ions (middle panel) 
and electrons (bottom panel) from R-3, in which o) ( |/Q ~0.8. 





Figure 5. Temporal evolution: (a) electrostatic wave amplitude, (b) Same as (a) but over an 
expanded horizontal scale, (c) Frequency spectrum of data plotted in (a) from R-3. 
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Figure 6. Temporal evolution of (a) the driving electric field amplitude E y (t) for the power-law 
spectra with the exponent a=-l. (b) the wave amplitude, (c) perpendicular temperatures of 
heavy and light ions and (d) frequency spectrum of the data plotted in (b) from R-4. 




Figure 7. Temporal evolution of the drifts of (a) light ions, (b)heavy ions (c) relative drifts 
between heavy and light ions due to their ExB drifts, (d) same as (c) but due to their polarization 
drifts, and ((e) drifts of electrons from R-4. The time period shown in these panels cover the 
large-amplitude burst of electrostatic noise in Figure 6a. 
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Figure 8. Same as Figure 6, but from R-5 in which the exponent in the power law spectrum 
a=-2. 


